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1. Introduction 

The standard method to calculate the scattering phase shift on the lattice is the Liischer's finite 
volume method It utilizes the shift of energy spectra in the finite periodic box to calculate the 
scattering phase. The formalism is restricted to the elastic scattering region. One attempted to 
extend it above inelastic thresholds [||]. Although it can lead to a constraint imposed on the matrix 
elements of S -matrix, S -matrix elements cannot be obtained separately. 

We propose a different method to extend it above the inelastic threshold. This is fulfilled 
by extending the procedure recently proposed by HAL-QCD Collaboration ffl]. We focus on 
the asymptotic behaviors of the equal-time Nambu-Bethe-Salpeter (NBS) wave functions near the 
spatial infinity, where the information of T-matrix is contained. We construct a coupled channel 
effective Schrodinger equation by introducing an energy-independent interaction kernel, so that 
it can generate all the NBS wave functions simultaneously. Once such an interaction kernel is 
constructed in the finite volume, S-matrix is obtained by solving the effective Schrodinger equation 
in the infinite volume. (The resulting wave functions are nothing but the NBS wave functions, 
which contain the information of T-matrix in their long distance part.) In this paper, we use NA- 
NL coupled system to formalize the procedure. 

Contents are organized as follows. In Sect. ^[ we consider the asymptotic long distance behav- 
ior of the NBS wave functions for NA-NT, coupled system. We will see that they contain T-matrix 
in their long distance part. In Sect. [|, we extend an effective Schrodinger equation for the coupled 
channel system by introducing an energy-independent interaction kernel. The extension is per- 
formed so that it can generate all the NBS wave functions simultaneously as solutions. In Sect. |], 
we derive the factorization formula, which plays a key role in constructing an effective Schrodinger 
equation in a coupled channel system in Sect. ||. 

2. Asymptotic behaviors of equal-time Nambu-Bethe-Salpeter (NBS) wave functions 
for a coupled channel system 

We consider NA-NL coupled system with ~ 940 MeV, ~ 1 1 15 MeV, and ~ 1 190 
MeV. For notational simplicity, we treat them as bosons. We are interested in the energy region 

mjy J r m\<E <niN + m\ + m z , (2.1) 

which is a combined region of (i) elastic NA region: + hia < E < + and (ii) NA-NL 
coupling region above the AT threshold: + < E < + m\ + m % . 

We consider two incoming states at the same energy E as \N(p)A(—p),in) and \N(q)L(—q),in). 
p and q denote the asymptotic momenta for NA and NL systems, respectively. They are related to 
E as 

E = ^m 2 +p 2 + ^m 2 A + p 2 = ^m 2 N + q 2 + ^m 2 +q 2 . (2.2) 
We define the NBS wave functions for an incoming state \N(p)A(—p),in) as 

WA#A(x h xr,E) = Z N 1/2 Z A l/2 (0\T[N(x l )A(x 2 )]\N(p)A(-p),in) (2.3) 
WmMxi,x 2 ;E) = Z N ll2 Z 1L {l2 {0\T[N{x l )Z{x 2 )}\N(p)A(-p),in), (2.4) 



2 



A coupled channel extension to the Liischer 's finite volume method 



Noriyoshi ISHII 



and the NBS wave functions for an incoming state \N '(q)L(-q) ,in) as 

Wa-M*UX2;E) = Z N l/2 Z A l/2 {0\T[N(x l )A(x 2 )]\N(q)U-q),in) (2.5) 
Ynwl(x u x 2 ;E) = Z N ll2 Z z ll \0\T[N{x Y )Z{x 2 )}\N{q)Z{-q),in), (2.6) 

where N(x), A(x) and £(x) denote local composite interpolating fields for nucleon, A and E 
baryons, respectively. Z#, Za and Zx denote the normalization factors involved in the limit N(x) — > 

1/2 1 /2 1 /2 

Zn N ut(x), A(x) ->■ Z A 7 A ouf (x), and E(x) ->■ Z E ' E 0[tf (x), respectively, as x -> +°°. By using the 
reduction formula, these NBS wave functions are related to S -matrix as 

(N(p\ )A(p' 2 ), out \N(p)A(-p) , in) (2.7) 
= disc. + J d 4 Xl d 4 x 2 e ip 'i Xi (Di + m 2 N ) e ip >- {U 2 + m\) Yna-na(xi ,x 2 ;E) 
(N(q[ )L(q' 2 ), out \N{p)A(-p) , in) (2.8) 
J d 4 x x d 4 x 2 e iq '^ (Di + m 2 N ) e iq '^ (D 2 + m|) Ynz-,na(xi,x 2 ;E), 

and 

(Nip', )A(p' 2 ) , out \N(q)L(-q) , in) (2.9) 

j d 4 x x d 4 x 2 e ip '' Xl (□ i + m 2 N ) e ip '^ (U 2 + m 2 A ) }jf NAiNL (xi ,x 2 ; E) 
(N{q\ )L(q' 2 ) , out \N(q)Z(-q) , in) (2. 10) 

= disc. + J d 4 Xl d 4 x 2 e^ Xl (n 1 +mj r )e i ^ Xl (□ 2 + m|) V^e^e(xi,x 2 ;£), 

where "disc." stands for disconnected terms, and = d 2 /dtf — V? denotes d'Alembert operator. 
These relations lead us to the following asymptotic behaviors of equal-time restrictions of NBS 
wave functions near spatial infinity (|x — y \ — > oo) as 

Wn/^nk (x-y;E) = lim \j/ NA -jfA (*, x ,y,yo = 0;E) (2.11) 

2 e ipr 

~ e' pr + - - - ^NAm (s) + • • • 

E N (p)+E A (p) pr 

Vnz-na{x-Y,E) = lim WNZ-NA(x,x ,y,y = 0;E) (2.12) 

9 e iqr 

^ny.-na{s) 1 , 



and 



E N (q)+Ez(q) ' qr 



WNK\Ni.{x-y;E) = ]imw NA . m: ( X , Xo ,y,y = 0;E) (2.13) 

2 e ipr 

^ -. . r ^NA-M-( s ) 1""- 

E N (p)+E A (p) pr 

YNL-Nz(x-f,E) = lim Wm:-M.(x,xo,y,yo = 0;E) (2.14) 

a: ^+0 

2 e iqr 
^ e' qr + — -, r~r~£r~7 \ $nwz(s) — + ••■, 
E N {q)+E z (q) qr 
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where 2? denotes the on-shell T-matrix, i.e., (f,out\i,in) = 1 + i(2n) 8 • (Pf — Pi)^f.i, E^(p) 



m 2 N + p 2 ,E A (p) = Jm^ + jP^^p) = yjml + p 2 , wdr = x-y. 
Derivations of these asymptotic behaviors of NBS wave functions are similar to single-channel 
case [Q>|5[ |6|]. Here, we give a brief example of ^na-,na^J\ p) case. 

YNA-MnE) =Z N l/2 Z A 1/2 (0\N(f)A(0)\N(p)A(-p),in) (2.15) 
= Z N 1/2 Z A l/2 f 3^* (0\N(r) \N(k)){N(k) | A(0) \N(p) A(-p) , in) +/(r) , 

where I(r) is referred to as the inelastic contribution, which corresponds to contributions not asso- 
ciated with single nucleon intermediate state. Since 1(7) does not propagate long distance, we will 
neglect this term below. Because of Eq. (2J_), we have an equality 

(N(k l )\A(0)\N( Pl )A(p 2 ),in) (2.16) 
= (0\A(0)a Ntin (h ) \N{ Pl )A{p 2 ),in) 

+iZ N l/2 f d 4 Xl e ik ^ (Di+>4) (0\T[N(x 1 )A(0)]\N(pi)A(p 2 ),in) 

= zf 2 ENiPl) ( 2 ,)3 g 3 (fel _ pi) +z y 2 IMf^i f^M , 

-{k l -p l -p 2 ) 2 + m A -is 
which is used to replace the second factor in the integrand of Eq. ( |2.15| ) as 

yw^^'+l e^Xw (2 ' 17) 



1 &(N(k)A(-k);N(p)A(-p)) 



' E A {k)-E N {k)+E N {p) +E A (p) E A (k)+E N (k)-E A (p)-E N {p)-ie 



Eq. ( |2. 1 ID is arrived at by performing the integration. The integration can be carried out by noticing 
that the propagating degrees of freedom come from the pole contribution, i.e., E A (k) +E^(k) ~ 
E A {p) +En(p), where the off-shell 2F can be replaced by the on-shell 2? ' . The derivations of 
Eq. (pi), Eq. ([2.1 3[) and Eq. (|2.14[) are similar. 



3. An extension of effective Schrodinger equation for a coupled channel system 

To construct an effective Schrodinger equation, we define the following K functions by multi- 
plying Helmholtz operators on the NBS wave functions as 

K NA - NA (x;E) = (A + f) Yna-,na(x;E) (3. 1) 

K N1L - NA (x;E) = (A + q 2 )\i/ Ni: -jv A (x;E), 



and 



K NKNZ (x;E) = (A + p 2 )yna-M^E) (3.2) 
Knz^z{x;E) = (A + q^xi/Nz-j^QciE). 
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These Helmholtz operators filter out the propagating degrees of freedom in the asymptotic forms 
of NBS wave functions. As a result, these K functions become localized objects, i.e., K is nonzero 
only within the range of interaction. We factorize these K functions as 



Kna-,na{x;E) Kn\-n^(x; E) 
Knz-na(x;E) Knz-js[z(x;E) 



(3.3) 



d 3 y 



WNA;NA(y,E) YNA;Nl.(y;E) 
YNZ;NA&E) YNZ;NI.(y>E) 



with an E -independent kernel U. (The proof of the factorization is given in Sect. Q) Since K is 
localized, U has to be localized as well. 

If such a factorization is possible, we have 



(A + p 2 ) yna-na(x;E) (A + p 2 ) y N Am& E . 
(A + q 2 ) Ynz-na(x;E) (A + q 2 ) Ynz-nz{x;E) 

U N A;NA(x,y) U NA - NI .(x,y) 

U N T.;NA(x,y) U NI .-Nz(x,y) 



(3.4) 



Vy 



YNA;NA&E) YNA-m&E) 



This implies that any linear combinations of NBS wave functions with arbitrary complex numbers 
a and /3 as 



Yna{x;E) = a\i/ NA - NA (x;E) + p\jf NA . NZ (x;E) 
¥nz{x;E) = aYNz-NA{x;E) + Pynz^l(x;E) 



(3.5) 



satisfy 



(A + p 2 ) Yna(x;E) 
(A + q 2 ) Ynz(x;E) 



d 3 y 



U N A;NA(x,y) U NA - N z(x,y) 
U N Z;NA(x,y) UNZ;Nz{x,y) 



Wa&E) 
¥nz(t,E) 



(3.6) 



which serves as an coupled channel extension of the effective Schrodinger equation. 
Several comments are in order. 



(i) Eq. (3i)) does not depend on a particular choice of boundary condition, since a and j8 
can be arbitrarily chosen. This is a desirable property, which makes it possible to determine the 
interaction kernel U in the finite volume. (The volume has to be sufficiently large compared to 
the range of the interaction.) Once U is determined in the finite volume, S -matrix can be obtained 



by solving Eq. (3.6) in the infinite volume. (Eq. (3.6) generates NBS wave functions as solutions, 
which contain the information of T-matrix in their long distance part.) 

(ii) The interaction kernel U is most generally a non-local object. As is described in Ref. [Q], 
a practical way to construct such a non-local interaction kernel U is to rely on the derivative expan- 
sion, which makes it possible to construct U order by order in the derivative expansion by using 
a finite number of NBS wave functions of varying energy E and angular momentum L. The most 
challenging step in fulfilling this procedure is the construction of NBS wave functions at excited 
energies with variational method. This turns out to be feasible. (See Ref. [^], where this procedure 
is applied to AA-NE-LL coupled system.) 
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(iii) In the non-relativistic limit, p 2 and q 2 are approximately related as q 2 / (2jUjve) — p 2 / (2/I.jva) + 
m A — mz, where fl^z and m^A denote the reduced masses for AT and NA systems, respectively. In 



this limit, Eq. ( p.6[ ) reduces to an eigenvalue problem used in Ref . [|2p : 



A 
2/ijvA 



+ UnA:NA 



UnA;NZ 




y NA {x;E) 




Wna 


{x;E)~ 


m z -m A - 2^ + Unz-,nz 






— ^NR 


Ynx 


x;E) _ 



(3.7) 

Unz;na/(2Hnz), Unz,nz 



where Una-j/a = Una-na/(2pIna), Una^vz = Una;nz/(2}x na ), Unz-,na 
= Unz-niI (2Hnz), and £nr = p 2 / {2/j, NA ). 

(iv) In constructing Eq. ([TfJ), we assumed the orthogonality of NBS wave functions. If the 
violation of orthogonality becomes severe, NBS wave functions should be orthogonalized before 
all the above procedures begin. The orthogonalization should be performed without affecting the 
long distance behaviors, where information of T-matrix is contained. 



4. Proof of the factorization of K function 

We assume that the following NBS wave functions are linearly independent 



WNA;NA{x',E 
YNZ;Na(x;E) 



YNA;NI.(x;E) 
YNZ;Nz(x;E) 



+m\<E<mN+mx+m ]t 



There exists a dual basis 

{ [VNA-,NA(x-,E),\jf NA . N z(x-,E)] , [YNZ;NA(x;E),\j/NZ;Nz(x;E)} \ 

I J mN-\-m^<E<mN+mh+m J[ 

which serves as 

Wna-na{y,E') WNA-m&E') 
Wnz-na{y,E') WNZ-Nzif> E ') ^ 

factorization can be verified in the following way: 



d 3 y 



a "left inverse" as an integration operator in the following sense: 

= (2k)8(E-E') 



YNA;NA(f,E) WNA-MZ&E) 
YNZ;NA(y;E) YNZ;NZ.(y> E ) 



dE' 
2k 



Now, the 

Kna-,na{x;E) Kna;nz{x;E) 
Knz-,na(x;E) Knz,nz(x;E^ 

Kna-na(x;E') Kn A ^nz{x;E') 
Knz-,na{x;E') Knz-nz(x;E') 

VNA;NA(3\E') WNA;Nz{f> E ') 

Vnz-ma&E') Wz-mz&E') 

Kna-na{x;E') Kn A; nz{x',E') 
Knz,na{x;E') Knz\nz{x;E') 

Wnama{y,E) Yna-nz(y e 
Ynz;na(t,E) Ynzmz&E 




¥na;na(t,E) Wna-mz&E) 
WNZtfA&E) Ynzmz&E) 

Wna-na{y,E') WNA-MziT'E') 
Wnz-na{y,E') Ynz;nz(yE') 



(4.1) 



(4.2) 



(4.3) 



(4.4) 
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We arrive at the factorization formula Eq. (|3.3[) by defining an interaction kernel U as 



U N A-NA(x,y) U NA - NI .(x,y) 



(4.5) 



dE' 
2n 



Kna-jva(x; E') K^A;NY.{x\E') 

Knl-na(x;E') Knz-,nz(x; E') 



WAm(y> E ') YNA;NI.(y;E') 

YNY,-N\{f,E') Wm.;NY.{f,E') 



Note that, due to the integration of E' , the interaction kernel U does not depend on the energy E. 



5. Summary 

We have proposed an extension of the Liischer's finite volume method above the inelastic 
threshold. We have considered NA-NL coupled system as an example. We have seen that equal- 
time Nambu-Befhe-Salpeter (NBS) wave functions contain the information of T-matrix in their 
long distance part. We have constructed an effective Schrodinger equation for a coupled channel 
system by introducing an energy-independent interaction kernel, so that it can generate all the NBS 
wave functions in the energy region + «2a < E < + mA + m n . Since this interaction kernel 
is localized, it can be constructed by using lattice QCD calculation in a finite volume. Once it is 
constructed in the finite volume, S-matrix can be obtained by solving the effective Schrodinger 
equation in the infinite volume. 

It is interesting to apply this method to the hyperon systems, where only a limited amount of 
experimental information is available. The first attempt to apply this method to AA-./VS-IX coupled 
system is reported in Ref. [0|, where the interaction kernel is constructed at the leading order of the 
derivative expansion by using NBS wave functions obtained by the variational method. 
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